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Abstract 

It is shown that the density of two initially independent condensates which 

are allowed to expand and overlap can show interferences as a function of time 

due to interparticle interaction. Using many-body theory, explicit expressions 

for the density are given which are exact in the weak interaction limit. General 

working equations are discussed which reproduce exactly the density in this 

limit. Illustrative examples are presented. 
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The investigation of interferences between particles is one of the most basic tools to learn 
on the nature of quantum gases. Interferences attracted much attention in particular in the 
case of Bose-Einstein condensates (BECs) both from the theoretical and experimental sides, 
see, e.g., [1-7]. In a popular set up studied, identical atoms are produced in two traps which 
we may call the left and right traps and which are separated by a barrier. By removing the 
traps and the barrier between them, the atoms expand freely and can overlap. In experiment, 
the photographs obtained show spectacular interference fringes [1,2]. 

The interference of two parts of a single coherent condensate is by now well understood, 
see, e.g., [4,8,9]. On the other hand, very little is known on the interference of two initially 
independent (i.e., fragmented) BECs, except for the case of non-interacting particles [3,5,6]. 
Fragmented BECs can be produced using a barrier between the two traps which is so high 
and broad that tunneling between them is negligible. 

In the available experiments, the atoms are prepared in a double-well trap potential and 
it is not generally proven whether the atoms form a coherent BEC, a fragmented BEC, or a 
combination thereof. Moreover, it is feasible nowadays to produce in the lab two spatially 
separate, initially independent BECs, see, e.g., [10], and this allows for experiments with 
definitely fragmented BECs. Apart from its importance as a fundamental problem, the 
solution of the problem of interference of two initially independent condensates is thus also 
of practical relevance. 

It is obvious from the literature that the impact of interaction on the interferences be- 
tween two initially independent BECs is not clarified. This is documented, e.g., in the 
review [7]. Can at all interactions between two initially independent BECs lead to visible 
interferences in the density? If yes, how do they look like? On what do they depend? These 
are the questions to which we provide clear and novel answers. 

In the scenario of two initially independent BECs the initial state of the many-body 
system before removing the traps reads 




(1) 
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where the b\ and b]^ are the usual creation operators for bosons in the left and right traps, 
respectively, which contain definite numbers A'^^: and Nji of atoms in them. After removing 
the traps, the state |^) is no longer an eigenstate of the system's Hamiltonian Hq and 
expands in space as a function of time. The time-dependent density, i.e., the expectation 
value of the density operator p{x) as a function of time becomes [7] 

p{x, t) = \p{x)\ = Nl \^l{x, t)f + Nr \^r{x, t)f , (2) 

where the $^^^(0;, t) are the single-atom states corresponding to 
bL,R{t) — Q^p{iHot)b^Rexp{—iHot). Obviously, the density is a sum of the individual den- 
sities of the two condensates and does not exhibit an interference term. 

We would like to draw attention to the fact that the literature result (2) has been obtained 
under the assumption that atoms belonging to the two different BECs do not interact with 
each other. In the following we demonstrate that in the presence of interaction, the density 
p{x, t) does show an interference term. This finding has many consequences. In particular, 
the corresponding interference structures remain after the statistical averaging over many 
experimental runs. Of course, as p{x, t) changes in time, the average must be carried out at 
the same value of t. 

The density operator p{x) — ^\x)^{x), where ^{x) is the usual field operator, can be 
expressed in any complete basis set of one-particle functions. Using (box normalized) plane 

waves, one has 

p(x)^v-'j:e'^'-''Kl,a,, (3) 

k,k' 

which is a suitable choice for our scenario. Here, v is the volume and is the destruction 
operator of a free boson with momentum k having the usual commutation relations [a^, a|./] = 
6kk'- In this representation, the Hamiltonian of interacting identical bosons of mass m, after 
release of the traps, is given by 

k 

^ = ^ '^fel+fc2,fe3+fc4«li«l2^fe3«fe4- (4) 

ki,k2,k3,k4 
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For the ease of presentation, the widely used contact interaction V{x,x') — XaS{x — x'), 
where Aq is proportional to the s-wave scattering length, is used [7,11]. Of course, any other 
interparticle interaction can be used as well. As usual, Hq describes the motion of the free 
atoms. 

We may now proceed to evaluate the time-dependent density p{x, t) of the interacting 
particles. First we transform to the interaction picture and write pi{x,t) = e*^°*p(a;)e~*^°* 
and all other quantities in this picture are extremely easy to evaluate. In particular, ak{t) = 
e*^°*ajte~*^°* = e~*2T?r*afe. To proceed we make use of the identity 

p{x,t)^{^i{t)\pi{x,t)\^i{t)) (5a) 

where \^i{t)) — e*^°*e~*-'^* |^(0)). This expression for the density can be systematically 
evaluated using the text-book expansion 

|^,(^)) = {l - ^ / Vj{h)dh + (-^)' Vi{h)dh J^' Vi{h)dt2 + . . .} 1^(0)) . (5b) 
Inserting this expression into the expression (5a) for p{x, t) gives 

p{x,t)=p\x,t)+z f' dti{^\[Vi{ti),pi{x,t)]\^) 
Jo 

+ fdh f dt2{-^\[Vi{t2)AVi{h)rpi{xM\'^) + --- (6) 

Jq Jo 

which is a systematic expansion of p{x,t) in terms of the interparticle interaction. |^) = 
1^(0)) is the initial state in Eq. (1), which is an eigenstate of the system in the presence of 
the trap potentials, and p^{x,i) = \pi{x,t) \ ^P) is the free density evolving without the 
impact of this interaction. With (6) and the commutator relations for the boson operators 
Ofc and a|./, one obtains an explicit expansion of p{x,t). With growing order, the terms 
of this expansion contain higher products of destruction and creation operators, a,l._^ak2, 
4i«l2«fe3«fc4> and so on. 

How can one now evaluate the desired density p{x,t) in view that the initial state |^) 
in Eq. (1) is expressed in terms of other operators and not in terms of a|.? To this end 
we remind that the field operator ^ can be expressed in any complete set of functions, 
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plane waves or others. Since bi and obviously commute, the respective functions $l and 
$R may be seen as the first two members in the expansion of the field operator: ^(x) = 
^L{x)bL + + "Rest". It immediately follows that 

ak = Mk)bL + Mk)bR + " Rest{k) (7) 

where (p{k) is the Fourier transform of ^[x) and the "Rest" contains all other destruction 
operators commuting with 6^ and hu. Luckily, all the latter operators have a vanishing 
impact when applied to the initial state |^). Consequently, if we order all the creation oper- 
ators to the left and all the annihilation operators to the right when taking the expectation 
value of p{x, t) with the many-body state |^) in Eq. (1), we can throw the "Rest" out in the 
following. Inserting the result (7) into the above mentioned higher products of destruction 
and annihilation operators in the expansion of p{x,t), one readily obtains 

p{x, t) = Pll{x, t) + Prr{x, t) + Plr{x, t), (8) 

where pll and prr are the densities of the expanding separated BECs as if the two BECs 
do not communicate, and p^R is the change of the density due to the interaction between 
them. The terms contributing to pll (prr) contain only hiihR) and h\{h^j^ operators, e.g., 
bibibLbL, and those contributing to plr contain only mixed products, e.g., b^ibRbibR- 

Note that the calculation of p{x, t) for interacting bosons amounts to solving a full many- 
body problem. Wishing to arrive at an exact analytical result, we concentrate first on weak 
interparticle interactions. 

With V in Eq. (4) and p{x) in Eq. (3) and the very simple appearance of ak{t) given 
above, we evaluated explicitly the commutator [Vi{ti), pi{x, t)] and with it the leading term 
in Eq. (6). The calculation is somewhat lengthy, but very straightforward. Using the simple 
"trick" (7) we obtain the exact results up to first order 

Plr = -4:XaNLNR Im{$i(x, t)ALR{x, t) + t)ARL{x, t)} , 

Pll = Nl \^L{x,t)f - 2A,iVi(iVi - l)lm{<^L{x,t)ALL{x,t)} . (9a) 
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Here ^L,Rix,t) is the freely expanding ^l,r{x) and the amphtude Alr reads 

V ' , , , ^0 



(9b) 



To obtain puR just interchange L and R in (9a). The quantities Arl, All and Arr are 
obtained analogously from (9b). The integration over ti in Eq. (9b) can, of course, be 
performed explicitly, but then the summations over the three momenta are more cumbersome 
to carry out. Let us briefly discuss the result (9). Clearly, the interference term plr vanishes 
for t — > 0. Furthermore, pLR{x,t) vanishes as expected if the atoms do not interact with 
each other (Aa — > 0). The interference term pLR{x,t) is enhanced by the product NlNr of 
the numbers of atoms in the two initial BECs. 

The above results make clear that the interaction between the particles gives rise to an 
interference term in the density of two initially independent BECs of identical bosons. Before 
presenting a numerical example we go one step further and pose the question whether we 
can formulate a mean-field theory which reproduces exactly the exact small result (9a). 
Such a theory would open the door for real applications. The standard mean-field leads to 
the well-known and widely used Gross-Pitaevskii equation [7,11]. For coherent states this 
equation gives exact results for small A^. Clearly, it is inapplicable to fragmented states 
(1). For fragmented states a more general multi-orbital mean-field theory has been recently 
derived [12]. In the present scenario two orbitals are involved and the respective time- 
dependent mean-field (TDMF(2)) takes on the appearance (for the general derivation of 
TDMF, see [13]): 

i'^L ^V[h + Xa{NL - 1) \^l\^ + 2XaNR \^r\^] ^Jl, 

i^R ^V[h + Xa{NR - 1) l^i^l' + 2A„7V^ ^r (10) 

where the initial conditions are ipL,Rix, t — 0) — ^L,Rix)- h is the usual one-particle Hamil- 
tonian (in our scenario just the kinetic energy operator) and V = 1— \ipL — \'4'r){'4'r\ 
is a projector which ensures orthonormalization of the orbitals ipL and ipR [13]. 
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We now prove that the TDMF(2) in (10) exactly reproduces the exact many-body small 



Xa result (9a). In TDMF(2) the density can be expresses by p{x,t) = NL\iljL{x,t)f + 



Nr lipnix, t)f. Since Xa in (9a) is taken to be small, we may write ipi = t) + Si/jl{x, t) 

and analogously for ipR. Inserting into the latter expression for the density and comparing 
with (9a), we immediately identify Si^l: 

St/JL = -iXa{NL - 1)AIl - i2XaNRAl^. (11) 

To obtain St/jR just interchange L and R. Of course, we still have to show that t/^L = ^l+Si/jl 
with dipL from Eq. (11) indeed fulfills the TDMF equations, i.e., is the solution of (10). 
Taking the derivative of + ^i^L with respect to time and using i^i — h^i, leads to 

iijL = h^L + Xa{NL - l)Al^ + 2XaNRAlj, . (12a) 

A^]^ and A^j^ can be deduced from Eq. (9b). The time t appears there twice. The deriva- 
tive with respect to the upper limit of the integral leaves us with a separable triple sum 
of terms like t'"^/^ e*''i'^e*''?/(^"'VL(^i) = ^L{x,t), which leads to and 
respectively. The derivative with respect to t in the exponential function gives an expres- 
sion which is identical to —hAlj^ and —hA}^j^, respectively (remember that h — — 2^^)- 
Collecting all terms yields 

= hijjL + [Xa{NL - 1) + 2XaNR\^Rf] $i (12b) 

which is, as the orthonormalization is of second order in A^, identical with the TDMF(2) in 
Eq. (10) to first order in the interaction strength Xa- 

In the following we present two illustrative numerical examples. In the first the inter- 
action is very weak and the initial single- atom functions ^l,r{x) are normalized Gaus- 
sians located at ixo, i.e., ^l,r{x) — (:^)^/^ exp{— (x ^ xo)^/a'^}. This is a realistic 
choice for harmonic traps. The Fourier transforms of these functions simply are </?l,r(^) = 
a^/^/(27r)^/'^exp(— /c^a^/4) exp(=i=iA;a;o). They evolve in time as ipL,R{k,t) = e~'^2^*ipL,R{k) 
and their back-transforms give ^L^R{x,t) which are normalized Gaussians expanding in 



space as a function of time and can be found in many elementary text books. For the 
ease of presentation our example is in one dimension. After transferring the summations in 
Eq. (9b) to integrals over the momenta, all three integrals over k^, k2 and ki can be carried 
out analytically using that e''^^'^^^^^^ dk — sfiilq for any complex quantities as long 
as Reg^ > 0. We have checked that this condition is fulfilled for all three integrals. The 
final result is of the form Aj^r{x^ t) = /q dtif{ti)e^^^'^°'^^'^\ where /(ti) and g{x, Xq, ti, t) are 
simple but lengthy algebraic expressions. 

Fig. 1 shows the interference term of the density pLR{x,t) for very weak atom-atom 
interaction. Without interaction pLR{x,t) = and the density is simply NL\^L{x,t)f + 
NR\^R{x,t)f. The interaction leads to plr for which we use expression (9). To simplify 
the discussion, we put m — a/2 — h — 1 and express x in units of a/2, t in units of (a/2)^m, 
and Aa in units of (^7^)3^ • Furthermore, Nl = Nr = 500, xq = 6 and Xa = 2.5 • 10"''. At 
t — one has piRix, 0) = and then it starts to grow as ~ at very short times and much 
faster later on. Plr{x-i t) exhibits an oscillatory behavior which changes as time proceeds. 

In the following we apply the TDMF theory. As a first step we compute PLR{x,t) for 
the case of weak interaction discussed above. The results are also depicted in Fig. 1 and 
seen to coincide with those obtained using the analytic expression (9). Next, we enlarge the 
interaction strength Aa. For coherent states the time-dependent Gross-Pitaevskii equation, 
which is exact in the weak interaction limit, has been demonstrated in many cases to be 
applicable for intermediate and stronger interactions [7,11]. Similarly, there is reason to 
expect that for fragmented states the TDMF theory, which has been proven above to be 
exact in the weak interaction limit, is applicable well beyond this limit. We mention that 
TDMF(l) is nothing but the time-dependent Gross-Pitaevskii equation. 

For the very weak interaction discussed above we have chosen as initial conditions Gaus- 
sians located at ±Xo. Now, as the interaction is increased to A^ = 0.1, we choose the 
respective solutions of the stationary Gross-Pitaevskii equation at this Xa as initial condi- 
tions to account for the interaction when the harmonic traps centered at ±Xo are released. 
In Fig. 2 the density p{x,t) computed using the TDMF(2) equations (10) is shown as a 
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function of time. As seen in the figure, at t — the density consists of two separated distri- 
butions centered at ixo- The traps are removed at this time and the distributions start to 
broaden and to overlap. At about t = 2.4 one begins to see impact of the interference term 
in the density which becomes strongly pronounced as time proceeds. 

We conclude that the density of two initially independent condensates which are allowed 
to overlap can show interference effects in the presence of interparticle interaction. The 
physics of so called fragmented states, like the state in Eq. (1), is generally very different 
from that of coherent states [14,15]. Coherent states of condensates have been extensively 
studied, mostly in the framework of the Gross-Pitaevskii equation [7,11]. A BEC in a 
coherent state can exhibit interference fringes even in the absence of interaction [7-9,16,17]. 
Take, for instance, the coherent state = {N\)-^/^{b^)^ \vac) with 6^ = (ftt. ^ 6^/^. 



This immediately leads to p'^"^{x,t) = ^ \^L{x,t) + ^R{x,t)\^ in the absence of interaction 
between the atoms. In analogy to Eq. (7) we can determine the interference term p'lr{x, t) — 
A^Re($2^R). For the expanding Gaussians discussed above, the oscillatory part of p^°^ is 
simply given by cos[K{t)x\ with K{t) = SxQ{t/m) / {a* + At^ /w?). This interference term is 
qualitatively different from that arising due to the interaction between the particles. Another 
important difference between pf^{x,t) and pLRix,t) worth mentioning is that the former 
depends on the relative phase between and ^r{x), while the latter does not depend 

on this phase. 

Whether in an experiment the initial state is coherent or fragmented depends on the 
experimental conditions. It is beyond the scope of this work to argue whether or not the 
initial state in the currently available experiments on interference is fragmented. It is also 
not our intention to take side in the ongoing debate on whether these experiments detect 
the density or higher-order correlation functions, although we tend to share the opinion 
of some researchers see, e.g., [4,8,9,16], that the density is measured. What we can state, 
is that if one measures the density of two freely expanding initially independent BECs, 
it will only show interferences in the presence of interaction. This leads to the following 
proposal for an experiment which makes use of the fact that nowadays one can vary the 



strength of the interaction between the atoms [18,19]. Two measurements are necessary. 
If the measurement with interaction shows interferences which disappear upon measuring 
with the interaction turned off, then (a) the initial state was a fragmented state and (b) the 
interaction is responsible for the interferences. 

The theory presented here is easily extendable to any kind of interparticle interaction. 
It is also easily extendable to the case where one does not let the two BECs expand freely 
by removing the traps completely. One may, e.g., remove only the barrier and let the BECs 
expand in the new global trap. Since the interference structures depend on the interaction, 
a wealth of effects can be expected by varying the interaction, the form of the individual 
traps and of the numbers Nl and of the particles. 
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FIG. 1. (Color online) The interference term p^pi for very weak interaction strength 
(Aa = 2.5 • 10^'^) at two values of t. The solid curves are computed using the analytic result 
(9). The dots are computed using the TDMF(2) equations (10). The quantities shown are dimen- 
sionless. 
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FIG. 2. (Color online) The density p{x,t) of two condensates of 500 atoms each for Aa = 0.1 
as a function of time computed with TDMF(2) (black) compared to the density piL + Prr of two 
BECs which do not interact with each other, each computed with the Gross-Pitaevskii equation 
(red). The quantities shown are dimensionless. For more details see text. 
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